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ABSTRACT 


Multispectral scanner data in twelve spectral channels, 
in the wavelength range 0.46 to 11.7 vm, acquired in July 
1971 for three flight lines, were analysed by applying 
automatic pattern recognition techniques. These twelve 
spectral channels were divided into four wavelength groups 
(111, W2, W3 and W4), each consisting of three wavelength 
groups — with respect to their estimated probability of 
correct classification <P ) — in discriminating agricultural 
cover types The same analysis was also done for the data 
acquired in August, to investigate the effect of time on 
these results. The effect Of deletion of each of the 
wavelength groups on P , in the subsets of one to nine 
channels, is given. Values of P for all possible 
combinations of wavelength groups, ^n the subsets of one to 
eleven channels are also given. 


1. INTRODUCTION 

Multispectral scanner (MSS) data were analysed in subsets » ' one to twelve 
spectral channels, in the wavelength range 0.46 to 11.7 ym, for selected 
flightlines of the 1971 Corn Blight Watch Experiment^. These t-J"^ve spectral 
channels were divided into four wavelength groups (Wl, W2, W3 ai i W4), each 
consisting of three spectral channels (Table I). The purpose of “his study was 
to determine the statistical separability of multispectral measurements from 
agricultural cover types for evaluation of these wavelength groups- The 
agricultural cover types selected were: corn, soybeans, green forage (hay & 
pasture), and forest., In particular, the objectives of the study were: (1) To 
study the effect of deletion of all possible combinations of the four 
wavelength groups •• Wl, W2, W3 and W4, on the statistical separability and 
corresponding estimated probability of correct classification <P ) of the 
agricul,tural cover types. (2) To develop a criterion for evaluation of a 
combination of wavelength groups, based on the estimation of the probability of 
correct classification obtained by using this combination of wavelength groups 
in discrimination of agricultural cover types. (3) To investigate the effect of 
time on these results. 

The literature review of the statistical separability of agricultural 
cover types was done by Kumar and Silva (1977)^. In addition, they analysed the 
multispectral scanner data in wavelength range 0.46 to 11.7 ym for three 
flightlines. They found that in the subsets of one to six spectral channels, 
the combination of wavelength regions (where V, N, M and T denote the visible, 
near infrared, middle infrared and thermal infrared wavelength regions, 
respectively): V, VM, VNM, VNMT, VVNMT, VVNMMT, respectively, were found to be 
the best choices for getting good overall statistical separability of the 
agricultural cover types for the data acquired on July 16 as well as August 12. 





An «ffort was made to explain these z*esults on the basis of spectral properties 
of agricultural cover types. The overall statistical separability of the 
agricultural cover types was found xo be greater for the data of August 12 than 
the data of July 16. Kumar^ (1977) did a further analysis of similar nature, to 
evaluate explicitly each spectral channel, each wavelength region. and all 
possible combinations of wavelength regions for statistical separability, in 
terms of estimated probability of correct classification for agricultural cover 
types. Deletion of the channel 7 (0.61 to 0.70 ym) reduced by about two 
percent. The deletion of each of the other channels caused no reduction, or 
less reduction, in the values of P^, as compared to this channel. The deletion 
of the spectral channels constituting the visible wavelength region caused more 
reduction in P^* as compared to the spectral channels constituting any of the 
other wavelength -regions . The deletron of the spectral channels of near 
infrared wavelength region caused relatively small changes in the values of 
In the subsets of one to six spectral channels, the combination of wavelength 
regions V, VM, VKT, VHMT, VVNMT, VVNMMT; and T, NT, VN7 or V>rT, VNMT, VVNMT or 
VMMMT, VVNMKT were found to be the best choices for the data of middle of July 
and middle of August respectively. 

Since the spectral channels were divided into four wavelength regions, the 
results obtained could be interpreted on. a physical basis; on the other hand, 
each wavelength region (V, N, M, T) is not represented equally, in the sense 
that each wavelength region does not have the same number of spectral channels. 
The author does not have data of a multi spectral scanner where each wavelength 
region has equal number of channels. However, the author felt a definate need 
for doing an analysis of similar nature by dividing the available twelve 
spectral channels (Table D into four wavelength groups — Wl, W2, W3 and W4, 
each wavelength group consisting of three channels, so that each wavelength 
group is .equally represented in the sense that it has an equal number of 
channels. In addition, this analysis is more complete in the following aspects: 
(1) The effect of deletion of all possible combinations of the four wavelength 
groups on the estimated value of Pg is given. (2) In addition to the average 
values of Pq, maximum as well as minimum values ox for all possible 
combinations of wavelength groups in the subsets of one to twelve spectral 
channels are give . 


METHOD OF ANALYSIS 

Multispectra: scanner data in twelve spectral channels in the wavelength 
range 0.46 to 11.' ym, collectec with an optical -mechanical scanner at 
altitudes of 914 r 2133 meters. -.3000 to 7000 feet) over Western Indiana were 
analysed by applying automatic t -^ ttern recognition techniques. The wavelength 
bands of these twelve spectral c tannels are given in Table I. The data of three 
selected flightlines, acquired i^ July of 1971, were analysed. Each of these 
three flightlines had fair or good amount of each of the four agricultural cover 
types: corn, soybeans, green forage and forest. These three flightlines were 
selected carefully so that these combined could be considered to be 
representative of the four agricultural cover types in the Western Indiana. 

Black and white photography and gray scale printouts of the spectral 
channels of the flightlines were used to aid in locating the boundaries of the 
fields on the Digital Display. Sufficient number of fields of each agricultural 
cover type were selected carefully so that they could be assumed to be 
representative of the flight line. 

Using the same three flightlines and twelve spectral channels, an 
i dentical analysis was performed on the data acquired in August of 1971, to 
study the effect of time on the statistical separability of agricultural cover 
types. The multispectral scanner data was acquired on both dates (July and 
August) between 10.30 a.m. and noon time (local solar time). In addition, 
these data were of good quality and free from problems like lack cf sufficient 
ground observations, excessive cloud cover, etc. The analysis was done for the 
data acquired in July and August, because corn and soybeans have reached their 
maximum vegetative growth by these times, and one month of time is sufficient 
for significant changes to occur in the spectral properties of agricultural 
cover types. The author wanted to avoid the analysis of data from late 
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September afterwards, because soybeans are harvested in Scpteaber-October, The 

f^uthor tried to keep all the variables other than time imiiform in the two 
.iJuIy and August) sets of data as far as possible* For exasiple, an effort was 
{-aade to select about the sasie field boundwies for the two sets of data. A 
!^tal of sore than fields taken irosi three flight lines were analysed* 

Xach field was treated as an Independent unit and the fields of the saine 
. 4 gr icultural cover type were put in the same class. The statistics algorithm 
was used to compute the mean vector and covariance matrix (mean and standard 
-deviation) of the classes* Histograms of the agricultural classes defined above 
^dre used to check unimodality of the statistical distributions in individual 
-ohannels. The classes were redefined to ^iminate distinct multiple modes. 
Divergence is defined for any two density functions. . In the case of normal 
variables with unequal covariance matrices, divergence in 4t spectral channels 
Cl, C 2 » •••» is given in terms of mean vectors and covariance matrices of 
the classes^ 

A modified form of the divergence D^, referred to as ••transformed 
-divergence**, has a behavior** more like the probability of correct 
rmlassification than the divergence, D* 

JHy * 2 {1 - exp(-D/8)> (1) 

Transf ormed divergence has been used throughout this study . 

Although divergence only provides a measure of the distance between two 
class densities, its use has been extended to the sulticlass case by taking the 
.average over all pairs'^. Let denote the divergence between classes i and j 

*of a certain flightline, then the average • divergence over all class pairs of 
tour classes (each agricultural cover was treated as a separate class) is given 

■'^AVG “ g [^®T12 * °T13 * ®T14 * %23 '* %24 * %J4^ 


-Let-Drjjjjj s-minimu* of 


%12’ ^113^ ^23’ ^24 ♦ ^34 


(3) 


Let superscri ts 1 and 2 w. th the symbol •*I>r** denote the values of 
transformed diverg. nee for the cata acquired in middle of July and middle of 
-‘■August respectively. Let ^MIN3 values of (see 

-eq. (3)) in first, second and third xul^tline ^ respectively, for the'data 
‘Acquired in middle of July. 

^AVG * t^AVGl * ^AVG2 ^VGs]] 


®TMIN * ®TMIH2’ ®TM1N3^ 


i « 1, 2 (4) 

, i * 1, 2 (S) 


^Assuming each agricultural class has a multivariate gaussian distribution, 
"^e fcatiire selection processor was used to find I>rAVG 
^11 possible combinations of one to twelve spectral channels^ out of the 
available^ twelve spectral channels* 
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subset of r spectral channels 
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Mximized over all possible subsets of r spectral channels out of the available 
twelve spectral channels. From the values of the average transformed 
divergence 9 classification accuracy can be reasonably predicted from the 
results of Swain et al. (1973)^. 

Table I gives the wavelength interval and the corresponding wavelength 
group of each of the twelve spectral channels. Tables II and III give the 
effect of deletion of each of the four wavelength groups and all possible 
combinations of the four wavelength groups respectively on .MAX’ 

terms of the corresponding estimated probability of correct classification 
(Fc>* 


To fulfill one of the main objectives of the study — evaluation of all 
possible combinations of wavelength groups — the following criterion is 
proposed. 

Each of 12 available channels of the multispectral scanner can be placed 
in one of the four wavelength groups Wl, W2, W3 and W4 (Table I). Thus, any 
combination of spectral channels can be called as the corresponding combination 
of the wavelength groups. For example, channel combination 1 4 7 10 is 

called "combination of wavelength groups W1 W2 W3 W4". For a given 
combination of wavelength groups, for example, W1 W2 W3 W4, and the 

corresponding value of were calculated using the curve of Swain and King^, 
for all possible combinations' of four spectral channels that constitute the 
combination W1 W2 W3 W4. The minimum, maximum and the mean of these values 
of Pq were calculated for all possible combinations of wavelength groups, in 
the subsets of one to twelve spectral channels, and is shown in Table, IV, for 
the data of middle of July as well as middle of August. 


3. KESULTS AND DISCUSSION 

Table II shows that deletion of each of the wavelength groups causes the 
following maximum reductions in P- for the data of. July: W1 (0.5, subset of 
five channels), W2 (0.4, subset of four channels), W3 (1.0, subset of one 
channel), W4 (2.3, subset of tw > channels). The corresponding values for 
reductions of I , for the data ci 'quired in August are: W1 (G.2, subs of nine 
channels) W2 (( 5.3 subset of n ne channels) W3 (0.5, subset ^f three 
channels), W4 (1.5, subset of o .e channel). 

Each wavelength region (visible, near infrared, middle infrared and 
thermal infrared) contains inde indent valuable information. Deletion of W4 
causes more reduction in P^ thu.. deletion of any of the other wavelength groups 
because it is the only group t'n-t contains all the spectral channels of two 
wavelength regions -- middle in: rared and thermal infrared. Deletion of W3 
causes more reduction in P^ thai' W1 or W2 because it has one spectral channel 
in the visible and two in the near infrared wavelength region, whereas W1 or W2 
has the spectral channels only in the visible wavelength region. 

Similarly Table III shows that, among the combinations of two wavelength 
groups, the deletion of W3 and W4 causes most reduction in because it has 
all the spectral channels of the near infrared, middle infrared and thermal 
infrared wavelength regions and one spectral channel of the visible wavelength 
region. Similarly, among the combinations of three wavelength groups, deletion 
of W2, W3 and W4 causes the most reduction* 

Table IV is quite useful, since it evaluates all possible combinations of 
wavelength groups in the subsets of one to eleven spectral channels, in terms 
of It shows that in the subset of one channel, W1 has the highest value 

of pQ* This is mainly because channel 7 (red) is an excellent channel for 
discriminating agricultural cover types. Besides other reasons pointed out by 
Kumar and Silva^ (1977), the separability of agricultural cover types in this 
channel was relatively high due to their low variance in this channel. In the 
subset of two to four spectral channels, the combination of wavelength groups 
W2 W4, W2 W3 W4 and W1 W3 (2W4), respectively, for the data of July; 
and W3 W4, W2 W3 W4 and W2 W3 (2W4), respectively, for the data of 
August, are foiind to be the best choices. This work has much application to 


UUJGiNAL fWGKII 

QB SQQ& sKttaaai 


f«atur« selection (i*e. selecting best subset of m spectrel channels, out of 
N available channels, of the existing multispectral scanners), and deciding 
wavelength bands of future satellites. 
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TABr: I. WAVELENGTH BANDS OF THE SPECTRAL CHAIiHELS 


Channel No. 

Wavelength Band 
(Micrometers) 

Wavelength Region 

Wavelength Group 

1 

0.46 - 

0.49 

visible 

1 

2 

0.48 - 

0.51 

visible 

A 

3 

0.50 - 

0.54 

visible 

1 

• 4 

0.52 - 

0.57 

visible 

2 

5 

0.54 - 

0.60 

visible 

2 

6 

0.58 - 

0.65 

visible 

2 

7 

0.61 - 

0.70 

visible 

3 

S 

0.72 - 

0.92 

near infrared 

3 

9 

1.00 - 

1.40 

near infrared 

3 

10 

1.50 - 

1.80 

middle infrared' 

4 

11 

2.00 - 

2.60 

middle infrared 

4 

12 

9.30 - 

11.70 

thermal infrared 

4 



: uTABLE II. EFFECT OF DELETION OF EACH WAVELENGTH GROUP ON THE 
PERCENTAGE OF CORRECT CLASSIFICATION 

' Number of 

chunels in Values of Probability of Correct Classification Estimated from 
the subset 



> 

0 


CM 

< 

^3 

A4 

®0 


»2 

®3 

®4 

1 

8^.^ 

84.3 

84.3 

83.3 

84.3 

85.4 

85.4 

85.4 

83. S 

85.4 

2 

89.2 

89.2 

89.2 

88.3 

86.9 

91.1 

91.1 

91.1 

90. C 

90.4 

3 

90.9 

90.9 

90.9 

90.2 

90.0 

94.3 

94.2 

94.3 

93.6 

93.5 

H 

92.6 

92.6 

92.3 

92.2 

90.6 

96.2 

96.2 

96.1 

95. t 

95.0 

s 

93.9 

.93.7 

93.7 

93.8 

91.8 

96.8 

96.8 

96.7 

96.7 

95.9 

6 

94.6 

94.2 

94,2 

94,4 

92.4 

97.3 

97.3 

97.1 

97.1 

96.2 

7 

9S.0 

94.6 

94.7 

94.7 

93.0 

97.5 

97.5 

97.3 

27.2 

96.4 

8 

95.4 

95. 0 

95.1 

94.9 

93.5 

97.7 

97.6 

97.5 

27.3 

96.5 

9 

95.7 

95.2 

95.4 

95.0 

93.7 

97.8 

97. C 

97. S 

97.3 

96.5 

10 

95.9 

np 

np 

np 

np 

97.6 

np 

np 

np 

np 

11 

96.0 

np 

np 

np 

np 

97.6 

np 

np 

np 

np 

12 

96.1 

np 

np 

np 

np 

97.9 

np 

np 

np 

np 


Note: This ;able gives i*ie ve.lues of percentage pr ^babij^ity of orrect 

classification estimated from the value** of (sc^ eq, ( 6 )), 

Swain and King C1973)^. Aq » Ai, A 2 * A 3 and A., denote *tne v ilues of 
when using all available channels, deleting >ectrai chan ‘ 'Is in the 
wavelength groups 1, 2, 3 and 4 respectively for the data acquired in 
the middle of July. Bq, I?2» B3 and 84 dem ote correspo* uing quantities 
as Aq* A]^, A 2 , A 3 and A 4 respectively for th data accuir* i in the middle 
of August. ”np** denotes that this channel co; bination was not possible. 
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TABLE III. EFFECT OF DELETIOK OF COMBINATION OF WAVELENGTH GROUPS 
ON THE PERCENTAGE OF CORRECT CLASSIFICATION 
SATA OF JULY 

(A) V«lu«t of PrabAbility of Correct Classificetion (Pc) Estimated from 

(see eq. (6) ) After Deletion of Combination of Wavelength Regions 


• 


V1V2 

W1W3 

V1W4 

W2W3 

W2W4 

W3W4 

W1W2W3 

W1W2W4 

W1W3W4 

W2W3W4 

1 

• 4,3 

84.3 

83.30 

84.3 

82.68 

84.3 

83.30 

82.68 

84.3 

83.30 

81.59 

2 

19.2 

89.2 

88.33 

86.38 

86.85 

85.79 

34.74 

86.12 

86.88 

86.06 

83.70 

3 

90.9 

90.9 

90.18 

86.96 

90.84 

89.62 

38.40 

88.56 

89.51 

86.45 

83.94 

4 

92.6 

96.26 

92.25 

90.50 

91.21 

90.87 

88.86 

np 

np 

np 

np 

S 

93.9 

92.89 

92.90 

91.35 

91.93 

91.64 

91.64 

np 

np 

np 

np 

6 

94.6 

93.89 

93.62 

92.29 

92.91 

92.68 

92.67 

np 

np 

np 

np 


SATA OF AUGUST 


(B) 

N 




Values of Probability of Correct Classification (P^) Estistated from 

(see eq. (6) ) After Deletion of Combination of Wavelength Regions 




V1W2 

V1W3 

W1W4 

W2W3 

W2W4 

1 

85.4 

85.4 

83.88 

85.4 

33.88 

84.05 

2 

9.1.1 

91.1 

90.95 

90.57 

90.01 

90.26 

3 

94.3 

91.64 

93.95 

93.55 

93.62 

92.33 

4 

96.2 

95.77 

95.96 

95.01 

95.16 

94.50 

5 

96.8 

96.44 

96.69 

96.26 

96.05 

95.88 

6 

97.3 

96.77 

96.77 

96.38 

96.18 

96.01 


W3W4 

W1W2W3 

W1W2W4 

W1W3W4 

W2W3W4 

93.51 

83.62 

83.88 

83.51 

83.51 

38.14 

90. 8S 

90.71 

38.93 

84.20 

90.53 

92.17 

92.10 

89.36 

85.47 

90. 7i 

np 

np 

np 

np 

92.05 

*-np 

np 

np 

np 

33.89 

np 

np 

np 

np 


Note: "np" denotes that deletion of combination of wavelength regions was not 
* possible. Wl, W2, W3 and W4 denote wavelength groups 2, 3 and 4 
respectively (see Table I). 

N « Number of channels in the subset 
Without Deletion 

c c 
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TABLE IV. LVALUA.'fOU OF COMBINATION OF WAVELENGTH GROUPS 
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